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HF  Prediction  Using  Adiabatic  Invariant  Theory 


I.  INTRODUCTION 

The  Propagation  Branch  of  the  Electromagnetic  Sciences  Division  of  HA  DC  is 
planning  to  conduct  an  experiment  which  will  explore  the  use  of  low-loss,  elevated, 
ionospheric  ducting  channels  in  order  to  achieve  super-long  range  propagation  at 
HF  (6-30  MHz;.  If  successfully  exploited,  these  propagation  channels  will  extend 
the  distances  reached  in  surveillance  and  communications  from  one  or  two  hop  to 
round-the-world  (RTW). 

The  purpose  of  the  experiment  is  to  collect  swept  frequency  ionograms  which 
will  provide  information  concerning  the  modes  of  the  super -long  range  propagation 
channels.  These  modes  will  be  combinations  of  the  three  idealized  modes 
shown  in  Figure  1.  The  first  mode  is  the  classical  ground -hop  mode  in  which  ravs 
bounce  back  and  forth  in  the  channel  between  the  earth  and  the  ionosphere.  The 
second  mode  is  the  chordal  mode  in  which  rays  are  refracted  by  the  ionosphere  at 
angles  close  enough  to  the  horizontal  such  that  thev  travel  along  chords 
connecting  points  in  ;  ie  the  ionospheric  shell  and  never  touch  the  earth.  The  third 
mode  is  the  trapped  mode  in  which  ravs  bounce  back  and  forth  in  one  ot  the  sparselv 
ionized  channels  between  the  relatively  denser  F,.  FI  and  F2  layers.  The  chordal 
and  trapped  modes  are  called  "earth -detached"  modes  because  their  ravs  never 
touch  the  earth.  Nevertheless,  thev  can  be  received  (or  launched)  from  the  ground 

(Received  for  publication  20  March  1084) 


in  the  vicinity  of  appropriately  tilted  ionospheric:  layers  that  refract  them  into 
(or  out  of)  the  ground-hop  mode. 


Figure  1.  Schematic 
Representation  of 
Three  Idealized  RTW 
Propagation  Modes 


2  GROUND-BASED  EXPERIMENTS 


Ground-based  long  range  and  RTW  propagation  experiments  have  provided 
evidence  that  ducting  takes  place.  ^  Figure  2  shows  a  RTW  ionogram  recorded  at 
2356  LT  during  one  experiment  which  involved  two  neighboring  ground  stations  at 
Ava  and  Verona,  New  York  (~  451 2  3  N  Lat.  )  that  were  only  65  km  apart.  ^  The 
137,5  msec  delay  time  of  the  leading  edge  indicates  that  the  signal  traveled  along 
the  longer  portion  of  the  connecting  North -South  great  circle.  In  this  successful 
pass  over  the  auroral  zone,  the  nighttime  auroral  trough  provided  a  natural 


1.  Bubenik,  D.  M.  ,  Fraser-Smith,  A.C.  ,  Sweeney.  L.  E.  ,  Jr.,  and  Villard,  O.  G.  ,  Jr. 

(197  1)  Observations  of  fine  structure  in  high-frequency  radio  signals  propagated 
around  the  world,  J.  Geophys.  Res.  761(4,):  1088 -  1092. 

2.  Elkins,  T.J.  ,  Toman,  K.  ,  and  Sales,  G.S.  ( 1980)  Theoretical  and  Experimental 

Studies  of  HF  Ducted  Propagation.  RADC-TR  -80 -"  SO,  AD  A098944. 
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i onus ;>!:<■  ri  c  till  i hat  guided  the  si  gna  1  up  an  i  me  r  ’  mo  :iti".)r;i  in  w ! i h  i :  nt  t ie : s. 
would  l;avt-  suffered  great  absorption.  Alth  High  i.'  • '  ■  > 1 1 .  I  -h  >:>  :■■(>  n*s  wen,  ■.  ert  ai  r.  i  >• 
present,  nevertheless  the  auriral  portions  he  path  n  vi.-.  ia\  !»••••.  -  '  .averse  i 
via  earth -detached  modes. 
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FREOUENCY  (MHil 

l-'i cure  2,  RTR  Ec-ho  of  KM'C  VV  Signal  Radiated 
Northward  l-’nm  Ava,  New  V-irk  and  Reeeived  at 
Verona,  New  York  About  137.5  msec  Later 

Some  features  of  this  ionograrti  are  wortn  noting.  The  entire  visible  :  rare  lie 
within  the  2.5  msee  time  interval  between  137.5  and  140  msec  delav.  The  broad 
vertieal  lines  which  run  from  top  to  bottom  below  13.5  Mil/  are  man-mane  inter¬ 
ference  and  noise.  Hevond  13.5  MHz.  where  the  interference  decreases,  the  train 
was  manually  changed  and  K  1AV  signals  were  received  out  to  25.  5  MHz.  At  the 
latter  frequency  the  gain  was  again  manually  changed  and  signals  were  seen  to  2!*  1 
The  splitting  of  the  trace  into  discrete  bands,  which  end  at  about  20.  23  and  211  MH 

•i 

is  tvpieal  of  RTVV  ionograms.  Tsedilina'’  has  theoretically  attribute  i  these  hands, 
whose  delav  times  increase  with  frequency,  to  ground-hop  and  chordal  modes.  On 
the  other  hand,  she  has  associated  other  band  structures  which  exhibit  decreased 
delav  times  with  increasing  frequence  with  trapped  modes.  However,  her  ,  onelu 
sions  are  baser!  upon  qualitative,  rather  than  quantitative,  matching  to  the  band 

3.  Tsedilina,  Ye.  Ye.  (1930)  Discrete  structure  oT  round-the-world  signals. 
Geomagnetism  and  Aeronomv.  20111:99-101. 
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structures.  The  availability  of  an  elevate;!,  orbitting  re 'fiver  in  the  HAIK'  experi¬ 
ment  would  remove  the  duet  ejection  losses  which  affect  ground-based  measure¬ 
ments  and  might  result  in  more  detailed  ionograms. 


.V  RADC  s  ORBITING  RECEIVED  EXPERIMENT 


KADC's  experiment  calls  for  signals  emitted  from  a  ground-based  source  to 
nter  one  of  several  possible  KTU  propagation  channels  either  via  a  naturally  tilted, 
ionospheric  refraction  gradient  (for  example,  the  solar  terminator'  or  via  a  nearby, 
artificially  produced,  ionospheric  scattering  region.  If  necessary,  the  latter  may 
be  created  by  means  of  the  heater''  shown  in  Figure  3,  The  receiver  will  orbit  the 
earth  every  90  min  with  perigee  altitude  within  the  ducting  region.  Nine  long-range 
ionograms  in  one  15-min  inverval  per  orbit  at  perigee  will  be  collected  and  telemetered 
to  earth  over  a  1-vr  period. 


PLASMA 

irregularities 


MAGNETIC  FIELD 


DUCTED 

MODE 


IONOSPEHERIC 

DUCT 


Figure  3.  Use  of  Heater-Induced,  F'ield-aligned  Plasma 
Irregularities  to  Deflect  Transmitted  Signals  Into  an 
Ionospheric  Duct 


The  success  of  the  experiment  depends  on  advance  knowledge  of  the  locations 
of  the  most  robust  ducting  channels.  Although  several  transmitter  sites  will  be 
available  during  the  experiment,  the  receiver  will  be  able  to  listen  to  only  one  of  them 
at  a  time.  Turn-around  time  for  ionogram  telemetry  and  processing  will  be  about 
three  days.  By  this  time  the  receiver  will  have  moved  through  many  orbits. 


Successful  scheduling  ol'  the  transmitters  will  require  quick  ami  reliable,  i  oiiqr:- 
tational  pre-assessment  of  the  relative  robustness  of  the  ducting  ch.mr.els  t<>  tin- 
receiver.  Since  conventional  ray  tracing  would  be  too  slow,  we  hue  e  turned  to  . 
technique  recently  developed  for  this  kind  of  problem  by  researchers  in  tin  soviet 

Union. 


4  ADIABATIC  INVARIANT  THEORY' 

-4  i  a 

The  idiabatic  invariant  technique  was  pioneered  by  Gurevich,  I'sedilina,  ’ 

g 

E'ishchuk  and  others  over  the  past  10  years.  It  is  an  outgrowth  of  earlier  work  hv 

YVovk.  '  The  theory  begins  with  the  magnetic  vector  potential  Air,  0  ,  o> 

'  l//(  rlV,  t£>,  O'j  which  is  inserted  into  Maxwell’s  equations.  Here,  r  is  the  radial 

distance  from  earth  center,  6  is  the  polar  angle  which  measures  distance  along  a 

great-circle  path,  o  is  the  off-path  azimuthal  angle,  4/  represents  the  radial  part 
^  — * 

of  A,  V ,  ^represents  the  angular  part  of  A,  and  j  is  an  arbitrary  unit  vector.  After 
separation  of  variable  and  linearization,  the  spatial  dependence  of  i*/i/.  =  r-  r  •  in  a 
spherically  symmetric  ionosphere  is  governed  by  the  one-dimensional  Schrocdingcr 

g 

equation  for  a  free-particle  in  a  potential  well: 

.2  2 

+  k  (E  -  U(z»]^  =  0  ,  (1) 

dz^ 

where 


E  =  f(f  +  l'/k2r^ 


U(z>  =  -  e  (z)  -  2z/r  , 
o 

€ (z>  =  1  -  fp(z)/f2  . 


rQ  is  the  earth  radius,  k  is  the  free-space  wave  number,  Mz)  is  the  ionospheric 
plasma  frequency  and  f  is  the  operating  frequency.  The  negative  of  the  potential  well 
function  U(z>  is  the  sum  of  the  magnetic  field-free,  collision-free  dielectric  con¬ 
stant  € ( z )  and  the  "spherical  distortion"  term  2z/r  ,  where  z  represents  altitude 
above  earth.  The  "energy  level"  E  represents  the  coupling  constant  to  the  differ¬ 
ential  equation  for  Y^. 


(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  37,  ) 


!•  lgure  4  snows  how  the  t  unction  t  i/.'  is  ■•oust  •  a  t«  ■■  i .  !'• .  tin-  Kit  is  .  ;  >!  ■ 

frequency  profile  over  the  point  '  1 ".  .4  s,  ii.  li  |;i  taken  from  the  I  I'S  K  »\(  \  1  • 
median  model  ionosphere  with  sunspot  number  4  0  in  .line  .t  11100  i  1'.  1'h,  .  j 
has  been  suppressed  for  clarity  .  To  the  right  is  'he  corresponding  I  :  >:•<  ■:'! !  • 

an  operating  frequency  of  20  Mil/.  Notice  that  the  latter  profile  :,lw.<vs  berm- 
C  =  -  1  on  the  ground  and  continues  along  the  tilted  line  -2/.  r  until  signife  .,;,t 
ionization  appears  at  an  altitude  of  ibout  80  km.  Thereafter,  this  line  ,  r  her 
is  modified  by  tne  LI  and  F  2- layer  ionizations.  If  tiiis  profile  is  rotaied  ui  d.  r 
as  in  Figure  4,  the  outlines  of  a  potential  well  is  immediate!-.  .pp  u-eni. 


PLASMA  FREQUENCY 
PROFILE 


POTENTIAL 

PROFILE 


figure  4.  1’lasma  Frequence  Profile  <>ver  the 
Point  (ll'.-V  S,  <1.8  F)  from  h  >Nt  Al*  With 
Sunspot  Number  40  in  dune  at  12  1  1',  end  the 
Corresponding  Potential  function  at  f-  20  Mil/ 


4,1  Propagation  Channel  Types 

4.  1.  1  GROL  ND-IIOP  CUANNKI. 


In  Figure  4,  the  vertical  portion  of  the  left-hand  wall,  which  extends  : 
L  =  -1  ta  P  <x  t  corresponds  to  the  perfectly  reflective  ground.  The  rial 
wall  stops  at  z,  where  L  attains  the  maximum  value  K,,.  Therefore,  tile 


or 
dial 
•va  v 


reversing  portion  of  the  hottomside  of  the  F- layer  lies  below  the  horizontal  iim 
U 

lying  above  U 


Fj.  because  this  line  represents  the  maximum  level  for  hound  solutions.  1. 


Fj  correspond  to  unbound  solutions  which  penetrate  through  tin 
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Ti  •ansmission  Channel,  Vol.  1,  Teelirnques  for  Analyzing  Ion 
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Figure  5.  Potential  W  ell  at  f=  20  MHz  for  the  Plasma 
Frequency  Profile  in  Figure  4  W  ith  )■ -Layer  Ground-hop 
Modes  in  Channel  I,  E- Layer  Chordal  Modes  in 
Channel  II,  FE-trapped  Modes  in  Channel  111,  and 
F-Layer  Chordal  Modes  in  Channel  IV.  The  E-  and 
F-Layer  maximum  electron  density  heights  are  h  E 

and  h  F,  respectively 
m 


ionosphere,  while  levels  lying  between  U  and  Ej.  and  U  =  - 1  correspond  to  bound 
solutions  which  bounce  back  and  forth  between  the  ground  and  the  ionosphere.  The 
latter  levels  belong  to  the  F-layer  ground-hop  channel  which  we  shall  call  Channel  1. 
Below  the  line  U  =  -1,  the  left-hand  wall  is  detached  from  the  ground.  Consequently 
we  shall  refer  to  the  corresponding  levels  as  earth-detached  levels. 

4.  1.  2  EA RTH- DETACHED  CHANNELS 

The  earth-detached  levels  are  divided  among  three  different  channels  due  to 
the  presence  of  the  secondary  maximum  F  =  Ej_.  at  z  .  From  U  =  - 1  down  to 
L  =  Ej_.,  F-layer  ionization  produces  the  right-hand  wall  while  the  upper  portion  of 
the  "pseudo  wall''  term  -2z/rQ,  which  contains  no  ionization,  accounts  for  the  left- 
hand  wall.  Therefore,  these  levels  correspond  to  the  F-layer  chordal  mode  channel 
which  we  will  call  Channel  IV.  Beneath  Channel  IV  lie  two  more  channels  both  of 


which  extend  downwards;  from  the  line  I  1.  .  i  h.mncl  II,  on  the  !•  ft,  is  the 
E -layer  chordal  mode  channel  which  lies  between  t  lower  portion  ■  >:  the  -2/  >• 
pseudo  wall  and  the  wave- reversing  portion  ol'  the  bottomside  of  the  K-l  ,\er.  l'o 
the  right,  between  the  topside  of  the  E-laver  and  the  low,-  portion  of  the  bottomside 
of  the  E-layer,  lies  Channel  111,  which  is  a  trapped  mode  channel.  If  an  I  1  bu-r 
were  present,  additional  channels  would  exist,  but  these  .'.ere  suppressed  in  order 
to  simplify  this  discussion. 

4.2  Snell’s  Law  in  the  Potential  Well 

Snell's  Law  in  a  spherical,  stratified  ionosphere  is  p  r  sin  i  -  const,  where 

P  =x/T  is  the  index  of  refraction  and  i  is  the  vertical  incidence  angle,  since 

t  ~  1  and  lit  =  -Jt  r-2  =  r  <  1 z  /  r  =  t  7t  i-l  t  2z  /  r  >  =  r  s/<  -  2/.  r  r  1  , 

_  v  ov  o  o  o  o  ov 

then  v - L  cos  a  s  const,  where  a  is  the  elevation  angle  complement;! ry  to  i  and  is 
measured  awav  from  the  horizontal  plane  passing  through  the  altitude  z  corres¬ 
ponding  to  E  =  L'(z).  Suppose  we  set  const  =  •J-X.  _  cos  ¥  at  the  reference 

altitude  z  ,  where  E  =  L'(z  >  represents  a  minimum  in  U.  Then  on  total 

c  r  e  j 

reflection,  a  -  0,  and  therefore  o  «  cos  y/ El  llT  gives  the  rav's  elevation 

angle  as  it  passes  through  z^ . 

For  rays  which  do  not  hit  the  ground,  a  -  0  occurs  at  two  different  altitudes 
which  define  the  upper  reflection  point  z^,  where  Lfz^)  =  E,  and  the  lower  reflection 
point  Zj,  where  U(z^)  =  E  also.  CJn  the  other  hand,  for  ground-hop  rays,  <v  =  0  is 

still  true  at  the  upper  reflection  point,  but  no  longer  on  the  ground.  Let  us  move 

the  reference  altitude  down  to  the  ground,  where  U  =  -1,  and  set  const  =  \/l  cos  , 
where  /3  is  the  ground  launch  or  arrival  angle.  Then,  a  =  0  yields  [i  -  cos 

We  can  use  a ,  z c  and  0  to  provide  reasonable  initial  guesses  for  rav-trace 
calculations.  More  significantly,  having  established  a  correspondence  between 
energy  levels  in  our  potential  well  and  rays  propagating  in  the  pertinent  ionospheric 
channels,  we  can  turn  to  the  problem  of  measuring  the  robustness  of  these  channels  - 
that  is,  counting  the  number  of  rays  (or  equivalently,  the  number  of  energy  levels' 
which  each  channel  can  accommodate. 

4.3  Number  of  Bound  Solutions  Inside  the  Potential  Well 

The  adiabatic  invariant  theory  enables  us  to  count  the  total  number  of  energy- 
levels  inside  the  potential  well  by  applying  the  WKB  solution  to  Eq.  ( 1'.  Outside  the 
well,  beyond  the  turning  points,  E  <  U(z)  and  this  solution  consists  of  damped  expo¬ 
nentials, 

4/(z)  =  1/2  C[iK(z)]  exp  [i  /  K(z'dz],  z  »  z  ,  and 

z 

u 
zc 

1/2  * 

exp  [i  /  K(z)dz],  z  «  z^. 


t^(z)  =  1  /  2  C'(iK(z)] 


(2) 


.‘.hen-  K'/.'-  K  v  i .  -  1  i/>,  l  .mil  l  i  iv  •  rlnt  I".  .  t.  •  n! 1  r  i-  i  r  .  '  :  ■  • , 

respoctivelv,  thi"  lower  and  upper  t'.irmnn  points  delim-d  1 : 1  - 1 1  ■  ■ : :  4.  !r-.rl«-  "1 

well,  between  the  turning  points,  tin -re  :uv  two  possible  k  mils  ■  i>- .  .:<•!  s..}:.u<.:i- 


,  .  -  l  ■* 

I  7.  I  t  '  |  K  /  >  J 


;  J  l\l  'll/  -  ~  4  | 


-  J  •) 

v/iz*  =  t":Kl/.>]  “  cos  |  /  Kizklz  -  jj  .’4],  z  ^  «.«.  /.  -v  ^  . 

'■( 


These  match  the  outside  solutions  onlv  if  k  1  =  tn  +  1  ■'  2 ' 

n 


1  -  /  v  E  -Uz>dz,  E  =  E  \  E.. 

n  ‘  n  o  n  \ 

ZC 


and  n  =  0,  1,  2,.  .  .  X  is  the  number  of  possible  nodes  in  U/  .  The  number  n  also 

n 

counts  the  number  of  eigenvalues  corresponding  to  bound  solutions.  We  see  that 
there  are  a  total  of  N+  1  eigenvalues  in  the  spherically  symmetric  ionosphere 
assumed  in  Eq.  ( 1 ). 

As  an  aside,  we  note  that  action  integrals  of  the  form  /  P  dz,  where 
P  =  V,2m(E  -  U)  =  phase  space  momentum  and 


phase-space  position. 


when  applied  to  periodic  systems  are  called  "adiabatic  invariants".  Hence,  the 
use  of  this  name  for  integrals  like  1^. 

4.4  Number  of  Bound  Solutions  in  Sub-Regions  of  the  Potential  Well 


In  Figure  6  we  show  I  vs  E  for  the  potential  well  of  Figure  5.  We  shall  drop  the 

2 

index  n  and  consider  I  and  E  to  be  continuous  because  1  s  10  to  10  km,  while 
X  =  2?r / k  s  10  m,  and  therefore,  kl  a  10^  to  10^  levels.  For  added  convenience,  we 
shall  ignore  the  multiplicative  constant  k  and  merely  let  1  represent  a  relative  number 
of  bound  solutions.  We  see  that  I  varies  continuously  from  the  upper  cutoff  (point  5, 
corresponding  to  E^.l  down  to  the  lower  cutoff  (point  3,  corresponding  to  E  l. 
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Figure  6.  Computed,  Relative  Adiabatic  Invariant  vs  Fnergy  Level  for  the 

Potential  Well  in  Figure  5 

From  Figure  5  we  know  that  this  curve  actually  spans  two  channels:  between  points 
5  and  4  it  lies  in  Channel  1,  and  between  points  4  and  3  it  lies  in  Channel  IV.  imme 
diately  below  point  3  there  is  a  gap  in  I  because  the  potential  well  splits  into  two 
subwells.  This  splitting  makes  I  a  double- valued  function  of  E:  the  curve  between 
points  1  and  2  lies  in  Channel  III,  while  the  curve  between  points  1'  and  2'  lies  in 
Channel  II.  The  gap  in  I  occurs  even  though  there  is  no  corresponding  discontinuity 
in  E.  We  will  discuss  the  implications  of  this  fact  in  Section  5.  For  the  moment, 
let  us  look  at  some  special  values  of  I  in  more  detail. 

4.4.1  ADIABATIC  INVARIANT  I 

The  number  of  bound  solutions  belonging  to  the  entire  well  is 


Where  E  in  Eq.  (4)  was  set  to  the  largest,  finite,  maximum  value  Ep.  and  the  inte¬ 
gration  runs  between  the  infinitely  high  ground  wall  at  z  =  0  and  the  altitude  z,_,  at 
which  U  =  Ep,  as  shown  in  Figure  5.  Recall  that  for  E  >  Ep,  solutions  can  escape 
over  the  right-hand  well  wall -that  is,  they  penetrate  through  the  ionosphere  and 
are  lost. 


4  4.2  \D1  \B\TU  l\\  \  HI  \\  r  !  , 


T!ii  -M-i  titui  nil  dr. tic  inv  .  r l  1 1’. !  i- 


i,  /  - 1  — • . 


where  I.  '  >>  set  u>  Hu-  \  ihif  -1,  .tI  im>-  ndegr  .tii>n  runs  between  tin-  e  "mind,  on 


a  :ni'h  I 


!,  tiio  lorn  >sj  4h- :i.  .Ptitude  ’■  .<t  w  liicli  tin’  multivalued  tunition  I  i.s 


..g:> in  <-c;no:  to  -1.  Hie  number  I  includes  bound  solutions  which  ;ut  :» 5 1  detached 
from  tilt-  forth,  since  I  •  -!  for  these  solutions.  It  olso  must  me  hide  tin-  grazing 

solution,  i'  I  t.  tn.it  hits  tile  ground  liuri/a  mta  I!  v.  If  tliis  grazing  solution  does 
not  exist  i  tint  is,  for  n  I  profile  with  L  «■  -  1  for  nil  /.  •  0»,  then  '.sill  be  defined 
equnl  to  zero.  This  does  not  mean  that  there  are  no  detached  solutions  in  this  case. 

I'o  till'  contrary,  tin-  only  possible  bound  solutions  would  indeed  be  detached  modes! 
However,  .«!!  of  these  solutions  would  already  have  been  included  in  I  ,  and  there¬ 
fore,  we  s'''  lj  d  to  indicate  this.  <  in  the  other  hand,  for  profiles  in  which  L 
varies  above  and  below  -t  tsuch  as  iti  Figure  5).  ^  will  ser\  e  r>  segragate  ground-hop 
solutions,  for  which  l  ■  -1  from  detached  mode  solutions,  for  which  l  •  -1. 

4.4.  .4  AD1A BATH.'  I.W  A  KI AN T  1 
The  third  adiabatic  invariant  is 


1  i  J  tlK  *  fhzldz  , 


where  If  was  set  to  the  secondary,  maximum  value  Ep,  and  the  integration  runs 
between  the  two  altitudes  and  z^.,  both  of  which  lie  beneath  the  E-layer  and  are 
intercepted  by  the  line  L  =  Ep.  Thus,  Ip  is  proportional  to  the  number  of  bound 
solutions  which  in  this  particular  case  are  chordal  modes  beneath  the  E -layer  in 
Channel  II  in  Figure  5. 

4.4.4  ADIABATIC  INV  A  RIANT  I,... 

1-  E 

The  fourth  adiabatic  invariant  is 


w  =  / 


Ep  -  U(z)dz  , 


where  Ep  is  the  same  secondary,  maximum  value  used  above  for  Ip,  and  the  inte¬ 
gration  runs  between  the  previously  mentioned  altitude  Zp  and  the  altitude  z^  where 
the  lower  portion  of  the  F-layer  is  intercepted  by  the  continuation  of  the  line  U  =  E  . 
Hence,  Ipp  is  proportional  to  the  number  of  trapped  mode  solutions  caught  between 
the  E  and  the  F -layers  in  Channel  III  in  Figure  5. 


»  ■  «-•  *  "A  •  %.  *  O  f  gU 


IT 


4.4.f>  ADI  A  HATH  1  \  \  A  K1 A  M’  I  (  ( ( . 

It  is  of  interest  to  note  that  the  same  encrg\  level  l.j  ’.‘.as  used  in  iletnuriL’  both 
lj.  and  I  .  Consequently,  the  latter  t  'A  <  >  may  be  summed  to  form  a  fifth  .uli.ib  itir 
invariant  Ij.^.  j.  -  1  »  1  which  represents  the  sum  of  the  bound  solutions  ..ssoeiateil 

with  the  presence  of  the  E-layer.  The  quantities  1.  ^  j_.  and  i(  ^  straddle  the  imp 
between  points  2  and  :i  in  figure  ti.  This  gap  shrinks  at  night  when  the  K-la\er  thins 
out  and  the  secondary  maximum  in  T  at  z^.  flattens.  In  general,  as  the  ionosphere 
becomes  depleted  of  electrons,  1  j 1  j . .  1|  |  >  and  lj  |  4  |.-  all  decrease  continuously. 
However,  as  noted  in  Section  4.  '.i.  2,  1  mav  abruptly  become  zero  if  the  ionosphere 
loses  enough  electrons  to  the  point  where  it  no  longer  supports  the  grazing  solution. 

4  5  Number  of  Bound  Solutions  in  Kach  Propagation  Channel 

The  adiabatic  invariants  Ij.  and  I  were  immediately  identified  with  the  number 

of  solutions  in  Channels  II  and  111  of  figure  a.  tin  the  other  hand,  it  is  possible  to 

count  the  number  of  ground-hop  solutions  in  Channel  1  by  subtracting  the  number  of 

earth-detached  solutions  1  from  the  total  number  of  solutions  lj  .  Likewise,  the 

number  of  f-layer  chordal  solutions  in  Channel  IV  may  be  counted  by  subtracting 

both  the  number  of  bl-laver  chordal  solutions  lj.  and  the  number  of  trapped  solutions 

I,.,,  from  I,.  Table  1  below  summarizes  the  "algebra'1  for  all  of  the  channels  at 
1-  E  1 

20  MHz  for  the  ionosphere  of  figure  4. 

I 


Table  1.  Propagation  Channel  Algebra  at  20  MHz  for  the 
Ionosphere  of  Figure  4 


Channel 

No. 

Channel 

Type 

Relative  Number 
of  solutions 

I 

F-layer  ground  hop 

V  -  *1 

II 

E-layer  chordal 

*E 

III 

E-layer/F-layer  trapped 

!fe 

IV 

F -layer  chordal 

!i  '  Ve+e 

At  this  point  we  the  important  result  that  we  have  estimated  the  number  of  bound 
solutions  within  each  channel  in  this  spherically  symmetric  ionosphere  without  ever 
having  traced  a  single  ray. 
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4. ft  Number  of  Bound  Solutions  in  Sub  Regions  ot  the 
Potential  Well  at  Lower  f  requencies 


The  foregoing  analysis  w  as  i  ..  rriod  out  for  dd  Mil/.,  a  fairl  .  hu  •.  : :  .. 

a  lower  frequency,  say  M  Mil/,  we  would  expect  to  find  i  deeper  p"te::u  A  ell 
taining  even  more  bound  solutions  than  at  20  Mil/.  Also,  we  would  e\;,,  ,  t  t  :  i  :i.  i 
some  ground-hop  solutions  below  the  K- layer  as  well  as  below  the  I  -lie:.  I'm - 
is  indeed  the  case,  as  shown  in  f  igure  7. 


Figure  7.  Potential  Well  at  f=  8  MHz  for  the  Plasma 
Frequency  Profile  in  Figure  4  With  F-Layer  Ground-hop 
Modes  in  Channel  I,  F-L.ayer  Chordal  Modes  in  Channel 
II,  FE-trapped  Modes  in  Channel  III,  and  E-Layer 
Ground-hop  Modes  in  Channel  V.  The  E-  and  F-layer 


maximum  electron  density  heights  are  h  E  and  h  F, 

O  ryj  J-VA  * 

respectively 
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4.6.1  PROPACA  HON  CHANNELS  A  l'  ,  ,\||!/ 

We  see  from  Figure  T  that  some  of  the  prop.it>.> tim.  •iiunn>ds  <iel'in>-<l  n  - 
tion  4.  1  had  to  be  revised.  For  example,  c  hannel  I  >1  -layer  ground-hop  solutions*, 
which  extended  at  20  MHz  from  E  down  to  I  ■  !.  n>>  .  only  goes  as  far  dov.  u  > s 
E ...  The  maximum  energy  level  limiting  Channel  II  <E-la\ or  chordal  solution-' 
now  lies  at  F  -1.  The  definition  of  Channel  111  <  trapped  solutions)  remains  un¬ 
changed,  but  Channel  I\  (F- layer  chordal  solutions)  has  completely  disappeared 
due  to  E-layer  screening.  We  also  created  a  new  Channel  \  to  describe  the  new 
E-layer  ground -hop  solutions  which  now  fill  the  region  from  E  down  to  C  -1. 

At  this  point,  there  seems  to  be  no  means  for  comparing  the  potential  wells  at  the 
two  frequencies  in  Figures  5  and  7. 

4.6.2  COMMONALITY  OF  THE  ADIABATIC  INVARIANTS 

In  spite  of  the  differences  in  well  topologv  at  the  two  frequencies,  the  formal 
definitions  given  in  Section  4.3  for  1^,,  1^.,  1^  and  1^  apply  to  both  wells.  The 
only  circumstantial  changes  involve  the  values  of  the  integration  limits  and  the  fixed 
constants  in  the  kernels.  Therefore,  we  can  still  use  the  four  adiabatic  invariants 
as  basis  elements  from  which  to  construct  a  new  channel  algebra  at  8  MHz  as  given 
in  Table  2. 


Table  2.  Propagation  Channel  Algebra  at  8  MHz  for  the 
Ionosphere  of  Figure  4 


Channel 

No. 

Channel 

Type 

Relative  Number 
of  Solutions 

I 

F-layer  ground-hop 

rF  "  XF  E+  E 

II 

E-layer  chordal 

III 

E-layer/F-layer  trapped 

Ve 

IV 

E- layer  ground  hop 

JE  '  *1 

We  could  go  on  to  examine  even  more  exotic  well  topologies  and  accumulate  (or 
eliminate)  additional  propagation  channels.  However,  we  would  find  that  the  four 
adiabatic  invariants  1^,,  1^.,  and  I  form  a  complete  set  of  elements  which  we 
can  use  in  linear  combinations  to  count  the  solutions  belonging  to  these  channels. 

We  shall  therefore  use  these  adiabatic  invariants  to  compare  different  ionospheres 
over  different  parts  of  the  globe  at  the  same  or  at  different  frequencies. 
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I  m-  m-egoing  analysis  applied  to  a  spherically  symmetric  ionosphere,  l  et  as 
extend  this  analysis  to  a  set  of  points  along  a  particular  path  around  the  world. 

1  igure  a  shows  the  variations  in  solar  zenith  angle  \,  layer  heights  and  critical 
frequencies  vs  angular  range  0  along  the  RTU  path  which  heads  due  east  from 
(42f  N,  1203  12)  at  0  -  0  and  passes  through  the  ionosphere  of  Figure  4  at 
0  240  degrees.  Figure  9  shows  the  sequential  values  of  the  overhead  adiabatic 

invariants  1^.,  I j.  p-+  ^  and  1]  which  were  computed  assuming  that  the  ionosphere  was 
locally  spherically  symmetric  at  each  point  0  along  the  path.  The  curves  were  made 
for  f=  20  MHz.  We  see  that  this  path  cannot  be  entirely  traversed  at  this  frequency 
via  ground-hop  modes,  because  we  find  path  segments  along  which  lf  vanishes. 

Recall  that  where  exists,  the  solutions  counted  by  1  are  divided  into  two  sets: 
ground-hop  solutions  with  lj  <  I  S  If.,  and  chordal  solutions  with  1  <  1  .  On  the 
other  hand,  where  I;  vanishes,  only  chordal  solutions  are  possible. 

The  lightly  hatched  regions  in  Figure  9  show  where  ground-hopping  is  permitted. 
Notice  that  it  is  forbidden  along  the  path  segment  between  (9=90  degrees  and 
0=  180  degrees.  Furthermore,  at  0  =  143  degrees  the  only  possible  propagation 
mode  consists  of  a  very  small  set  of  chordal  solutions  lying  in  an  extremely  shallow 
earth-detached  well.  This  shallow  well  controls  the  entire  RTW  propagation  channel 
for  this  path  at  20  Mflz  because  all  ground-hop  and  many  chordal  and  trapped  rays 
penetrate  through  it  and  are  lost.  As  we  can  see,  the  adiabatic  invariant  theory 
gives  us  at  least  a  qualitative  picture  of  the  survivable  propagation  modes.  The 
meaning  of  1(0)  can  be  extended  as  follows. 

The  fundamental  theorem  of  the  adiabatic  invariant  approximation  states  that 

the  bound  solution  0  )  associated  with  l^lz.0  )  =  const,  propagates  in  a  slowly 

varying  ionosphere  in  accordance  with  Snell's  Law;  conversely,  (//^(z ,0  )  undergoes 

a  transition  to  another  \p  (z,0  )  associated  with  another  1  (9)  =  const.  '  only  via  a 

m  m  J 

sudden,  scattering  process  which  violates  Snell's  Law.  A  formal  proof  of  this 
theorem  is  not  yet  available.  The  ionosphere  has  to  be  slowly  varying  because  we 
will  assume  that  it  is  locally  symmetric  about  the  radius  through  any  particular 
surface  coordinates  (0  , <2> )  for  the  spatial  extent  of  at  least  one  period  of  ray  oscilla¬ 
tion,  9  .  By  ’  period  of  oscillation"  we  mean  the  distance  interval  required  for  a  ray 
to  depart  and  return  to  the  same  layer.  Examples  are  the  great-circle  lengths  of 
either  one  hop  in  a  ground  -hop  channel,  one  complete  bounce  out  and  back  between 
the  two  layers  of  a  trapped  mode  channel,  or  one  ricochet  from  and  back  to  the 
single  layer  of  a  chordal  mode  channel. 
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Figure  8.  Variation  in  Solar- 
Zenith  Angle,  Critical  F requencies 
and  Maximum  Electron  Density 
Heights  (from  lONCAP  with 
Sunspot  Number  50  in  June  at  12  UT1 
Along  a  Great  Circle  Path  Heading 
Due  East  from  (42°  N,  120°  E'  and 
Passing  Through  ( 19.  5 0 S,  6. 8° E> 
at  240  Degrees  Angular  Range 
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Figure  9.  Adiabatic  Invariants  Ip,  Ij,  lpE  and  IFE+E  vs  Angular 
Range  at  f  =  20  MHz  Along  the  Great  Circle  Path  in  Figure  8 
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W  *.•  have  already  established  the  relationship  between  Snells'  l.aw  and  E  ,  and 

therefore,  1  ,  in  Section  4.  2  for  a  symmetric  ionosphere,  lienee,  in  the  symmetric 
n  - 

case,  we  have  a  unique  ray  which  traverses  the  altitude  /.  at  the  local  elevation 
angle  a  corresponding  to  1  .  We  shall  assume  that  this  is  still  approximately 
true  in  the  nearly  symmetric  case,  and  that  E^ifL,  z^(6),  and  z  (tf) 

readjust  in  such  a  way  as  to  leave  I  id  )  invariant. 

Therefore,  the  adiabatic  invariant  approach  allows  us  to  "rav -trace"0  simply 
by  superimposing  horizontal  lines  corresponding  to  constant  values  of  I (t9  1  upon 
Figure  9  as  shown  in  Figure  10.  By  moving  either  eastward  or  westward  along  one 
of  these  lines  we  can  immediately  determine  the  propagation  modes  of  the  corres¬ 
ponding  ray.  For  example,  the  segment  A^AAjAg,  along  which  1(6)  =  2.5,  des¬ 
cribes  a  nearly  horizontal  ray  that  ricochets  in  a  chordal  mode  near  the  bottom  of 
the  nighttime  F- layer  since  KtB  =2.5  corresponds  to  very  low-lving  levels  E (6)  and 
very  shallow  elevation  angles  om9).  In  daytime  along  AgA^,  this  ray  becomes 
trapped  between  the  E-  and  the  F- layers  since  2.  5  <  The  foregoing  is 

similarly  true  for  all  rays  within  the  band  I  (6)  i  2.  5. 


0  60  120  180  240  300  360 


ANGULAR  RANGE  8  (d*g) 

Figure  10.  The  Adiabatic  Invariants  Chart  ih  Figure  9  With  Raypaths 
AAjAgAgA,  BBjBg,  CCjCg.  DD ^DgDgD^DgDg  and  LL^Lg  (shown  as 

horizontal  dashed  lines)  Straddling  the  Longer  Great  Circle  Segment 
Between  the  Vertical  Dashed  Lines  at  9  =  90  and  240  Degrees 
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On  the  other  hand,  rays  with  !•(>>  •  2.  .*>  leave  the  ionosphere  somewhere  aheie 
this  l{T\\  path.  For  example,  ray  H  leaves  at  B.;,  either  by  directly  penetraline. 
through  the  1-  -layer  or  by  first  bouncing  off  the  ground  and  then  penetrating  some¬ 
what  farther  along  the  path.  This  uncertainty,  which  could  be  resolved  hv  .u  ta  ills 
doing  the  raytracing,  is  of  no  concern  as  long  as  we  know  that  this  ray  is  ultimately 
lost.  Indeed,  the  advantage  of  the  adiabatic  invariant  thoorv  lies  in  its  ability  to 
deal  with  the  question  of  existence  more  directly  than  the  constructive  approach  of 
ray- tracing. 

To  examine  a  more  practical  case,  suppose  we  are  located  on  the  ground  near 
9  =  240  degrees  and  we  wish  to  transmit  to  an  orbitting  receiver  near  9  -  t'O  degrees. 
We  first  draw  two  vertical  dashed  lines,  at  d  =  240  degrees  and  9  =  90  degrees  in 
Figure  10,  and  then,  begin  a  process  of  eliminating  rays  (in  other  words,  horizontal 
lines)  which  emanate  from  the  vertical  line  at  0  =240  degrees  until  we  find  an  unimpeded 
propagation  channel.  Since  we  wish  to  transmit  from  the  ground,  only  line  segment 
DC'Bl.K,  which  lies  in  the  ground-hop  channel,  is  of  any  use.  All  horizontal  lines 
pointing  westward  from  DCBI.K  penetrate  through  the  ionosphere  before  a  =  1 B 0 •  . 
Therefore,  we  need  to  consider  only  eastward  pointing  rays.  Of  the  latter,  only  the 
rays  lying  below  LL^l^  reach  0=90  degrees;  those  above  1,1,  L,  penetrate  through 
the  F-layer  before  9  =90  degrees.  Hence,  we  are  left  with  segment  lX'BL.  We 
might  also  consider  the  option  of  eliminating  the  rays  between  lines  BB^B^  and 
LLjLg,  because  they  involve  some  ground-hopping  between  BB(  and  LI^  that 
increases  the  likelihood  of  D-  and  E-laver  absorption. 

The  rays  between  lines  DD  ^  and  CCj  are  controversial  because  they  eventually 

fall  into  the  gap  (darkened  region)  between  lj.^.  and  1°  Section  4.4  it  was 

pointed  out  that  I  underwent  a  sudden  discontinuity  within  this  gap,  although  the 

level  E  did  not.  There  are  two  possible  points  of  view  regarding  this  situation.  The 
10  5 

liberal  view  ’  maintains  that  there  is  a  possibility  that  in  traversing  point  Dj  the 
ray  continues  to  obey  Snell's  Law  while  gradually  converting  from  chordal  to  trapped 
mode.  This  process  is  schematically  represented  by  the  vertical  line  It 

then  proceeds  in  trapped  mode  from  to  D^,  where  it  encounters  another  gap  in 
I  and  reverts  to  chordal  mode,  in  a  process  represented  by  D^D^.  Finally  the  rav 
continues  along  D^D^Dg  in  a  chordai  mode.  On  the  other  hand,  a  conservative  view 
would  hold  that  all  of  this  is  unlikely,  and  that  the  only  indisputably  open,  low 
absorption  channel  lies  between  lines  CCjCg  and  Perhaps  the  KADC  experi¬ 

ment  will  shed  light  on  this  issue. 

10.  Lobachevskiy,  L.  A..  Suares,  J.  ,  Tushentsova,  LA.  Fishchuk,  D.  I.  ,  and 
Tsedilina,  Ye.  Ye.  (1979)  Analysis  of  the  conditions  of  propagation  of 
decametric  radio  waves  on  a  long  path  by  the  adiabatic  method  on  a  fixed 
frequency.  I,  Geomagnetism  and  Aeronomy,  19(3):314— 319. 


Cv 

lU 

I" 


1  •• 


A  more  exhaustive  analysis  of  this  communication  channel  would  require  ua  to 
redraw  Figure  10  at  Frequencies  above  and  below  20  MHz  in  order  to  find  the 
optimum  frequency  which  would  make  this  channel  most  robust  — that  is,  the  fre¬ 
quency  which  yields  the  widest  band  of  acceptable  1  values.  Also,  we  should 
generate  additional  adiabatic  invariant  charts  for  other  times  of  day,  seasons, 
sunspot  numbers,  and  so  on.  Because  of  lack  of  space,  we  will  consider  only  one 
more  case  at  f  =  8  MHz. 

From  Figure  1  1  we  see  that  at  8  MHz  there  exists  an  earth-detached  channel, 
similar  to  the  one  at  20  MHz,  below  line  AA  A.  It  also  consists  partly  of  1  -layer 
chordal  modes  and  partly  of  trapped  modes.  However,  at  8  MHz  the  F- layer  is 


1-  L 


more  influential  and  the  trapped  modes  predominate.  Compare  the  longer  1 
curve  in  f  igure  11  with  the  shorter  curve  in  f  igure  9.  In  addition  to  this 
RTW  channel,  we  find  another,  which  consists  purely  of  ground-hop  rays,  between 
lines  CCjC.jCgC  and  BB^B.^B^B.  lie  might  consider  adopting  this  channel  for 
communicating,  either  eastward  or  westward,  between  9  =  240  degrees  and 
6=  90  degrees  (vertical  dashed  lines  in  f  igure  11)  via  ground-hopping.  The  westward 
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igure  11.  Adiabatic  Invariants  I(.,  Ij, 


'ff  and  Vf+F  vs  Angular 


Range  at  f=  8  MHz  Along  the  Great  Circle  Path  in  Figure  8,  With 
RTW  Haypaths  AA^A,  BB^B^B^B,  CCjCgCgC  Shown  as  Horizontal 

Dashed  Lines.  The  shorter  ravpaths  BB0,  CC^  and  DD  straddle 
the  shorter  great  circle  segment  between“the  vertical  dashed  lines 
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direction  would  be  preferable  because  it  is  shorter  .tin!  ao  .>  ;!  !  m  !u  h-  ulditi"n.,l 
rays  between  DIJ^  and  CV.j.  However,  li-  ,.nd  II  —  !:*>  •  *:•  ..b.-.irption  wo  in  : 
preclude  any  successful  g round-hopping.  Tnerefore,  we  see  that,  <•  •:  Alii.'., 
ground  to  elevated  receiver  communication  would  be  difficult.  Aa.ua,  it  is  important 
to  realize  that  we  have  arrived  at  these  conclusions  without  n  iving  done  .:n  .<■  tuil 
ray  -  tracing. 


6.  SIMULATED  IONOGRAMS 


If  we  prepare  additional  charts,  similar  to  Figures  !*  and  1  1,  for  other  fre¬ 
quencies  in  the  HF  band  and  superimpose  upon  them  horizontal  lines  of  constant 
I,  then  we  can  construct  a  simulated  ionogram  between  two  points,  using  a  procedure 
described  by  Fishchuk  and  Tsedilina.  6  For  each  frequency  we  determine  the  mini¬ 
mum  I  and  the  maximum  1  for  the  path  in  question.  We  consider  onlv  1  values 
min  max 

between  the  two  lines  1(d)  =  I  and  l(d>  =  I  .  For  each  intermediate  line  lit)),  we 

min  max 

compute  the  corresponding  curves  F(f)',  z^!d)  and  z^it<>  by  working  backwards  from 
Eq.  (4).  This  task  is  easy  to  do  it  we  can  interpolate  within  a  precomputed  table  of 
these  quantities.  We  then  insert  the  appropriate  values  into  Borisov's  1  formulas 
for  group  delay: 


z 


2  ,  zdz 

r0H  z(  y  E  -  U 


d  8 

yr-E 


(9) 


(10) 


and 


H  = 


dz 


/  - - 

ze 


(in 


where  t  is  the  group  delay  in  msec,  c  is  the  speed  of  light  in  km /msec,  rQ  is  earth 
radius  in  km,  dj.  is  the  path  end  point  in  radians,  z  is  the  ray  "mean  altitude"  in  km 
and  H  is  the  period  of  ray  oscillation  in  radians. 


11.  Borisov,  N.  D.  ,  Fishchuk,  D.  I. ,  and  Tsedilina,  Ye.  Ye.  (1978)  Frequency- 

dependence  of  the  group  lag  time  of  round-the-world  signals,  Radiophvsics 
and  Quantum  Electronics,  21(3):255-260. 
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At  this  point  we  have  mdisc nnunateh  counted  nil  of  tile  r  s  borne  bv  the 
ionosphere  between  9=0  and  9  =  9  ...  We  need  to  eliminate  those  rays  which  do  not 
emanate  from  the  source  and  those  that  do  not  arrive  at  the  receiver.  To  do  this, 
we  keep  track  of  the  number  of  oscillations  experienced  by  each  rav  via 

flf 

1>  =  /  dfl/t-i,  (  12) 

0 

and  we  accept  only  those  rays  for  which  P  is  an  integer.  In  this  way  the  simulated 
ionogram  will  represent  a  point-to-point  transmission. 

The  simulated,  lossless  ionogram  in  l-'igure  12  is  a  preliminary  result  and  will 
be  refined  in  future  work.  We  see  signatures  of  rays  which  undergo  from  14  to  33 
hops  KTW  lying  within  the  trailing  and  leading  edge  (dashed  curves)  for  ground-hop 
propagation.  These  results  are  similar  to  those  first  obtained  by  l-'ishchuk  and 

g 

Tsedelina.  In  addition  to  these  ground-hop  details,  we  have  also  computed  the 
trapped  and  chordal  mode  structures.  These  appear  as  two  bands  of  simulated 
signatures:  one,  immediately  below  140  msec  and  at  frequencies  less  than  12.  5  MHz, 
includes  rays  which  undergo  from  14  to  20  periods  of  oscillation  KTW;  while  the 
other,  immediately  above  137  msec,  includes  from  20  to  26  periods.  These  bands 

merge  into  one  band  which  lies  between  138  and  139  msec  from  12.  5  to  33  MHz. 

3 

These  results  do  not  match  those  of  Tsedelina  for  trapped  modes,  as  far  as  delay 
times  and  structure  are  concerned. 
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Figure  12.  Simulated,  Absorptionless,  Swept-frequency, 
RTW  Ionogram  from  an  Elevated  Sourch  Which  is  Radiating 
Due  East  from  {42°  N,  120°  E) 
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It  is  premature  to  compare  Figure  12  with  the  empirical  ionog ram  in  Figure  3. 
The  paths  and  ionospheric  conditions  were  not  the  same,  and  the  earth-detached 
mode  structure  as  shown  in  Figure  12  needs  more  work.  However,  we  can  find 
nothing  in  Figure  12  which  resembles  the  three  prominences  above  20  MHz  in 
Figure  3.  3  sedelina  qualitatively  attributed  these  prominences  to  surviving  ground  - 

hop  modes  but  her  quantitative  results  predicted  that  thev  should  lie  below  13  MHz. 


7.  ADIABATIC  INVARIANT  SYNOPTIC  MAPS 


Ionograms  are  suitable  for  predicting  operating  frequencies  over  the  path  be¬ 
tween  two  given  points.  For  the  purpose  of  assessing  propagation  robustness  between 
a  fixed  point  and  points  in  a  variety  of  azimuthal  directions,  a  synoptic  view  is 

preferable.  Here  is  where  the  real  power  of  the  adiabatic  invariant  technique  was 

12 

demonstrated  by  Tushentsova. 

Figures  13  and  14  are  synoptic  maps  of  f^E  and  f  F2,  respectively,  using  the 

IONCAP  model  ionosphere  for  June  with  sunspot  number  50  at  1200  UT,  in  a  format 

13 

similar  to  the  one  used  by  N  BS  to  show  critical  frequency  contours.  For  visual  im  pact, 
bands  of  adjacent  critical  frequencies  were  shaded  in  tones  of  gray,  instead  of  labelling 
each  contour  tine.  Also  in  this  format.  Figures  15  through  22  show  synoptic  maps  of  Ip, 

Ip,  Ip  and  Ipp  at  two  frequencies —8  and  20  MHz.  The  gray  tones  in  the  latter  figures 
indicate  bands  of  adjacent  I  values.  Some  unlabelled  contours  were  drawn  inside 
the  broader  shaded  bands  in  order  to  avoid  blurring  out  various  topological  details. 
The  dot-dash  curves  in  Figures  13  through  22  represent  the  solar  terminator,  while 
the  dashed  curves  denote  the  RTW  path  discussed  in  Section  5.  By  applying  the  over¬ 
lay  in  Figure  23  to  Figures  13  through  22,  variations  in  the  ionosphere  and  in  the 
related  propagation  modes  along  great-circle  paths  in  all  azimuthal  directions  that 
emanate  from  the  starting  point  at  (42°  N,  120°  E),  may  be  compared. 


12.  Tushentsova,  I.  A.  ,  Fishchuk,  D.  I.  ,  and  Tsedilina,  Ye.  Ye.  (1974)  Maps  of 

the  global  distribution  of  some  parameters  of  ionospheric  wave  ducts. 
Preprint  No.  5,  IZMIRAN. 

13.  Zacharisen,  D.  H.  World  maps  of  F2  critical  frequencies  and  maximum  usable 

frequency  factors,  NBS  Tech.  Note  2  (Apr  1959)  and  NBS  Tech.  Note  2/2 
(Oct  1960). 
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7  1  Synoptic  lp  Mapc 

The  maps  of  I  ,  (Figure  15  for  li  Mil/,  and  !  iguro  Hi  and  20  .Mil/)  slmv  tin- 
relative  robustness  of  111'  propagation  via  ground-hop  and  or  chord  >1  nodt-s  ur.  o!  .  me 
the  bottom  Side  of  the  global  H- layer.  The  contour  regions  in  tii.-.-o-  I  i:  .  .in¬ 

form  to  those  of  f  E  in  Figure  15.  To  further  distinguish  bct.-.-cn  g  round- ho;.  ,nd 
chordal  modes,  it  would  be  necessary  to  compare  these  map.-.  ith  tho-e  b-r  I  m 
Figures  I!'  and  20,  as  will  be  done  in  Section  7.5. 

7.2  Synoptic  Ip  Maps 

The  maps  of  Ip.  (Figure  17  for  li  MHz  and  Figure  18  for  20  Mll/i  sun.1,  the  rel  ,- 
five  robustness  of  total  HF  propagation  anywhere  over  the  globe,  irrespective  of 
the  mode  type  or  the  ionospheric  layer  involved.  The  I  contours  resemhle  those 
for  f  F 2  in  Figure  14;  the  most  robust  regions  occur  in  the  vicinity  of  the  Appleton 
anomalies,  where  the  1'2-layer  is  densest  and  highest;  while  the  least  robust  re¬ 
gions  occur  at  the  poles  and  in  the  nighttime  hemisphere.  The  effects  of  tin  auroral 
zones  are  not  visible  because  IONCAP  lacks  a  realistic  model  auroral  ionosphere. 

The  robustness  of  the  groundhop  and/or  chordal  modes  that  exclusively  involve  the 
F -layer,  (without  in  anyway  involving  the  E-layer)  may  be  measured  by  adding  the 
Ip,  and  Ip,p,  maps  together  and  then  subtracting  this  sum  from  the  corresponding 
Ip  map  at  the  same  frequency. 

7.3  Synoptic  Ij  Maps 

The  maps  of  Ip  (Figure  19  for  8  MHz  and  Figure  20  for  20  MHz)  show  where 
chordal  modes  exist  either  by  themselves  (Ip  =  0)  or  together  with  ground-hop  modes 
(Ip  7*  0).  While  ip  may  vary  continuously  within  the  gray  regions  (Ip  >  0),  the  transi¬ 
tion  from  gray  to  white  (Ip  =  0)  occurs  abruptly  at  points  where  the  ionosphere  has 
just  become  depleted  of  sufficient  electrons  to  refract  the  grazing  solutions  back  to 
earth,  as  was  discussed  in  Section  4.4.2. 

By  subtracting  Ip  (only  in  the  regions  where  Ip  4  0)  from  Ip.  for  the  corresponding 
geographical  areas,  it  is  possible  to  obtain  a  map  of  the  robustness  of  the  E-layer 
ground-hop  modes.  It  is  also  possible  to  obtain  a  similar  map  of  exclusively  F- layer 
ground-hop  modes,  by  first  subtracting  both  I£  and  I  from  Ip.,  and  then,  sub¬ 
tracting  Ip  (only  in  the  regions  where  Ip  4  0)  from  the  resulting  map.  Of  course, 
ground-hopping  is  forbidden  in  the  remaining  regions  where  Ip  =  0. 
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Figure  15.  Synoptic  Map  of  Adiabatic  Invariant  I..  at  f=  8  MHz  for  the 
IONCAP  Model  Ionosphere  With  Sunspot  Number  50  in  June  at  12  UT 
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Figure  16.  Synoptic  Map  of  Adiabatic  Invariant  Ijr  at  f=  20  MHz  for  the 
IONCAP  Model  Ionosphere  With  Sunspot  Number  50  in  June  at  12  UT 
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Figure  19.  Synoptic  Map  of  Adiabatic  Invariant  Ij  at  f=  8  MHz  for  the 
IONCAP  Model  Ionosphere  With  Sunspot  Number  50  in  June  at  12  UT 
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Figure  20.  Synoptic  Map  of  Adiabatic  Invariant  1^  at  f=  20  MHz  for  the 
IONCAP  Model  Ionosphere  With  Sunspot  Number  50  in  June  at  12  UT 
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7  4  Synoptic  lpj:  Maps 

The  maps  of  I  (Figure  21  for  [i  MHz  and  1  igure  22  for  20  MHz)  show  the 
robustness  of  trapped  mode  propagation  between  the  F-  and  1  -layers.  By  comparing 
Figure  21  with  Figure  15  for  f^E,  we  see  that  at  fi  MHz  this  mode  is  primarily 
controlled  by  the  E-layer;  while  comparing  f  igure  22  with  f  igures  14  and  15, 
indicates  that  at  20  MHz  both  the  E-  and  the  F-layers  are  influential.  Figure  22 
also  illustrates  how  the  RTW  path  discussed  in  Sections  5  and  fi  was  chosen  to  pass 
through  regions  that  were  more  robust  in  trapped  mode  solutions  at  20  MHz.  A 
North-South  path  through  the  60°  E  meridion  would  probably  have  been  as  good  or 
better,  but  was  not  used  because  it  would  pass  through  auroral  zones,  which  were 
not  realistically  modelled  in  lONFAP. 

8.  CONCLUSION 

We  have  shown  how  the  adiabatic  invariant  technique  was  developed  and  applied 
to  the  problem  of  describing  the  propagation  modes  along  great-circle  paths  (or 
path  segments)  in  order  to  qualitatively  assess  either  ground-to-ground,  ground- 
to-satellite  or  satellite-to-sateltite  communications.  We  have  demonstrated  how 
to  obtain  quantitative  results  in  the  form  of  simulated  ionograms.  We  have  also 
indicated  how  it  is  possible  to  obtain  a  better  understanding  of  global  propagation 
characteristics  via  synoptic  maps  of  1^.,  1^.,  and  1  .  We  have  discussed 
neither  how  to  incorporate  absorption  effects  into  ionogram  simulation  nor  how  to 
prepare  synoptic  maps  of  absorption.  This  will  be  undertaken  in  future  work. 
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Figure  21.  Synoptic  Map  of  Adiabatic  Invariant  1 
IONC'AP  Model  Ionosphere  With  Sunspot  Number 
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Figure  22.  Synoptic  Map  of  Adiabatic  Invariant  I__  at  f=  20  MHz  for  the 
IONCAP  Model  Ionosphere  With  Sunspot  Number  *  aO  in  June  at  12  UT 
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